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Abstract: We study the effect of an (9(A^^) density of heavy quarks in strongly 
coupled A/" = 4 SUSY Yang-Mills theory in the large N limit. This is achieved in the 
type IIB supergravity dual by introducing a uniformly smeared density of macro- 
scopic string sources stretching to the boundary of AdS^ x S^. The backreacted 
system exhibits a flow from an AdS^ "hedgehog" geometry to a scaling Lifshitz-like 
solution Lif 5 x with dynamical critical exponent z = 7, wherein the scaling symme- 
try is broken by a logarithmic running dilaton. We find an exact black brane solution 
within the scaling regime which describes the low temperature thermodynamics of 
the system. 



1. Introduction 



The AdS/CFT correspondence has been instrumental in unravelhng many 
facets of certain types of strongly interacting field theories in the large- limit. 
Within this framework, the behaviour of matter at finite density has been intensely 
studied with a view towards applications in condensed matter physics (see reviews p, 
^, 1^ and references therein). The physics of cold, dense QCD matter presents equally 
exciting challenges ||^, |^, |^ for holographic approaches involving flavour branes in 
both the probe limit @, |lO|, 11, 12| and incorporating their backreaction [^. 



The focus of this paper will be to study the effect of a sizeable, spatially homo- 
geneous number density of heavy quark "impurities" introduced as external sources 
into A/" = 4 supersymmetric Yang-Mills theory in the large- limit and at strong 't 
Hooft coupling. The number density of such heavy quarks will be chosen to scale as 
A^^, so that we may then study the backreaction of these on the gauge theory. One of 
the motivations for considering this setup is to obtain intuition for the finite density 
thermodynamics of holographic models with a large number of flavours (with Nf/N 
fixed in the large- A^ limit). 

According to the dictionary provided by gauge/string duality, a heavy quark 
transforming in the fundamental representation of the SU{N) gauge group corre- 
sponds to a macroscopic fundamental string stretching to the boundary of the dual 
geometry |jl^, |15[. The end-point of the string is anchored to the worldline of the 
quark in the boundary gauge theory. For a single quark, it suffices to treat the string 
as a probe in the dual geometry. It is natural to ask what happens when a large 
number of such quarks is introduced and, in particular, when this number scales 
as A^^ in the large- A^ limit, so that the quark impurities necessarily backreact on 
the gauge degrees of freedom. This is a difficult question to answer in general, but 
it is conceivable that a configuration preserving sufficient amount of global symme- 
tries may render the problem tractable within the dual gravity description at strong 
coupling. 

For a single heavy quark, or a straight timelike Wilson line, the corresponding 
string world-sheet is localized at a point on the transverse in AdS^ x S^, breaking 
the 5*0(6) global symmetry to 5*0(5). If we consider a large number of such string 
sources, we may distribute them uniformly along all transverse directions, thereby 
restoring both spatial translational invariance and the 5*0(6) global symmetry in 
the boundary theory. The energy density associated to such smeared distributions ^ 
scales as ~ nFi/(27ra') = npi \/X/2t[, where ripi is the number density of macroscopic 
strings or heavy quarks and A, the (large) 't Hooft coupling. The heavy quark 



^Recently, smeared string configurations have appeared in the study of finite density physics in 
a different system Mn. 
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distribution then backreacts on the degrees of freedom of the A/" = 4 theory when 

One might worry that a smearing of string sources along all transverse directions 
would lead to violation of Gauss' law applied to the sources. Indeed, the Chern- 
Simons couplings of IIB supergavity require that such smearing is consistent provided 
that a certain three-form flux is also switched on. Intuitively, the resulting three- 
form flux may be interpreted as a density of baryon vertices |T^ which soak up the 
charge corresponding to the string sources. 

Our primary goal in this work is to understand the zero/low temperature be- 
haviour of this highly symmetric smeared system as a simplified model for a finite 
density state in a theory with a known string dual. We further restrict attention 
to the A/" = 4 theory on M^'^. We postpone the study of the thermodynamics of 
the model, and specifically its finite volume phase structure, to future work [|T^. 
Since the quarks we consider are infinitely heavy, quark number cannot fluctuate 
and therefore all thermodynamic questions are restricted to the canonical ensemble. 
Related to this point is the absence of a global U{1)b baryon number symmetry in 
this system. However, there is a natural embedding of the above configuration into a 
holographic model with dynamical massive quarks and a U{1)b symmetry - namely 



the D3/D7 system with a large number of smeared flavour branes |T3|, [19 



The main result of this paper is that the backreaction of the heavy quark density 
on A/" = 4 theory triggers a flow to an infrared regime which exhibits approximate 
dynamical or Lifshitz scaling i.e. invariance under the scale transformations t — )■ 
t and x ^ ax. The dynamical critical exponent for the system under investigation 
has a fixed value, z = 7, independent of any free parameters ^ . Importantly, the 
Lifshitz scaling symmetry is only approximate, due to logarithmic running of the 
dilaton in the gravity dual. This behaviour closely resembles the scaling solutions 



found in ^ ^ ^ for charged dilaton black branes. As in the latter systems, 
the running of the dilaton renders a' corrections important deep in the infrared (IR). 
By electric-magnetic duality of the A/" = 4 theory, and S-duahty of type IIB string 
theory, the same statements apply to magnetic sources (replacing the strings by Dl- 
branes), except that now the dilaton becomes large in the deep infrared, necessitating 
the inclusion of string loop corrections. 

We also find an exact black brane solution in the Lifshitz-like scaling regime, 
which, by virtue of a small non-vanishing temperature, shields the system from the 
regime where the classical supergravity approximation breaks down. Taken at face 
value this solution has vanishing entropy density in the zero temperature limit, a 
physically appealing feature. However, clearly, in the strict zero temperature limit 



^After this paper was released, we were made aware by the authors of |2^] that this particular 
scaling solution has appeared in that paper (see Appendix C of pl|). 
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the corrections mentioned above will alter the conclusions somewhat (e.g. ^|). It 



is nevertheless extremely interesting to learn that the effect of a finite heavy quark 
density is qualitatively different to the effect of a finite R-charge density in the Af = 4 
theory. As is well known, the latter, for a specific choice of R-symmetry chemical 
potentials yields the Reissner- Nordstrom solution that has a large entropy density 
at zero temperature. It is worth noting that the emergence of Lifshitz-like scaling as a 
result of backreacting distributions of D-branes/baryons has also been seen generally 
in 



28 



It should be pointed out that the backreaction of heavy quarks or Wilson lines 
on A/" = 4 SYM has been studied in two different contexts. The first of these 
are the bubbling geometries dual to supersymmetric half-BPS Wilson lines in large 
representations (with ranks of order A^^) [p9| , pO| . The second context which is closely 



related to the discussion in this paper is the work of Headrick in [31 1 where the effect 



of smeared Polyakov loops on A/" = 4 SYM at finite temperature was investigated. 

This paper is organized as follows. In Section 2, we derive the reduced 5D action 
describing the smeared system of strings. The features of the ultraviolet behaviour of 
the resulting solutions are discussed and clarified in Section 3. Sections 4 and 5 are 
respectively devoted to the derivation of the IR scahng solution and the (numerical) 
construction of the flow from AdS^ x to the Lifshitz-like scaling regime. We discuss 
possible future directions in the final section. 



2. IIB plus strings 

As explained above, we take a finite heavy quark density in the field theory to cor- 
respond to a geometry sourced by a uniform distribution of fundamental strings, 
each stretching to the conformal boundary of AdS^ x S^. For simplicity the strings 
are chosen to be smeared uniformly along the internal 5*^ directions. This preserves 
the S0{6) global R-symmetry of A/" = 4 SYM, but breaks all supersymmetry. The 
coupled system of type IIB supergravity and a uniform density of macroscopic fun- 
damental strings is described by the action 

If 1 

S = Sub + Sfi H npi / B2 A dx^ A dx^ Adx^ Au^ — . (2.1) 

27ra' J TT'^ 

Here Spi is the Nambu-Goto action for the string sources, npi is the number density 
of heavy quarks (or strings) in the boundary theory. Go's the volume form on the unit 
five-sphere while B2 is the Neveu-Schwarz two-form potential which couples to the 
string world-sheet . 

We then look for a consistent solution to the supergravity equations which pre- 
serves 5*0(6) global symmetry, the spatial translational and rotational invariance. 
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and with the correct fluxes turned on to account for the presence of a non- vanishing 
fundamental string charge. The relevant field equations for the three-form and five- 
form fluxes i^s, F3 and ^5 = ^5-!- |i?2 A F3 — ^H^ A C2 in Einstein frame, are: 

d*F, = H^AFs, d{eKF^) = -g^HsAF,, (2.2) 

gs d (e"* * H3) = -gl F5 A F3 + — A dx^ A dx^ A dx^ . 

zira TT'^ 

The asymptotic value of the dilaton e'^ is set to gg the string coupling, while 

is Newton's constant in ten dimensions, with IGtcGn = {2Tiy gla'^. The field 

equations are then solved by 

(7,F5 = 4(l + *)a;5, i^3 = 0, (2.3) 

F3 = — p= dx A dx A dx . 

We have used the standard AdS/CFT dictionary, namely that the 't Hooft cou- 
pling A = RXds/^'^ ~ (^TTfy'siV), and we have conveniently rescaled the AdS radius 
-^Ads to unity. The introduction of macroscopic fundamental strings stretching from 
the boundary of AdSr, to the origin, necessarily requires a non-zero three-form field 
strength. 

A non-vanishing flux for F3 would suggest the presence of D5-brane sources. 
Although there were no explicit D5-branes in the setup to begin with, we could 
interpret the source of the three-form flux as a distribution of D5-branes wrapped on 
with a number density, ~ ^ . This is consistent with the well known fact that 
N quarks must be bound into a gauge-invariant baryon operator which corresponds 
to a wrapped D5-brane with strings attached [|I7|. Below we want to obtain the 
backreacted geometry wherein the fundamental string distribution and associated 
baryons can be replaced with smooth backgrounds accompanied by corresponding 
fluxes, and no additional D-branes. Since the fluxes are completely determined as 
above in terms of npi and A^, the only 'active' fields in the background are the dilaton 
and the metric components. 

2.1 Metric Ansatz 

Given the symmetries of the problem, the metric and the dilaton have the form 

ds"" = g„ir) dr^ + gu{r) dt" + e^-^^") dx^ + e'^'^'^^ d^l , (2.4) 
= e^g:\ 
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where we have defined the shifted dilaton (f) which vanishes at infinity. In this back- 
ground, the Nambu-Goto action for the smeared distribution of strings is, 

^Fi = J d'x J drV^^.e"'/'. (2.5) 

This is accompanied by a term quadratic in npi, arising from the energy density 
in F3. We will find it useful to work with the 5D effective action after appropriate 
rescaling of the lOD metric: 

9ap = e^'^'^/^g^fs, = 0,1,...4. (2.6) 



(is(5) = gap dx^dx^ = gu dt^ + cjrr dr"^ + e^" dx^ , = a + ^rj 



This ansatz along with eq.( |2.3|) , when substituted into the type JIB supergravity 
action, yields the reduced 5D system: 



S = -— d''x^g[R-fr^''r-\<P"r+ (2.7) 



The first of the potential energy terms is the scalar curvature of the . The con- 
tributions from the F3 and F5 fluxes, and the smeared Nambu-Goto action for the 
strings can be packaged neatly as a perfect square ^ . The bulk action depends on 
the quark number density via a single parameter Q defined as 

Q,VA.^ = VA.^. (2.8) 

Q is a density and therefore a dimensionful parameter. With regard to its dependence 
on the number of colours A^, we will treat it as a number of order one in the large 
limit, so that the quark density is 0{N'^), or equivalently, the baryon number density 
ns is 0{N). The system we are studying and the form of the action above is closely 
related to that of which did not need to include a non-zero F3 ^. 



3. UV AdS-hedgehog solution 

The backreaction of a uniform distribution of string sources on pure gravity (in four 
dimensions) was studied in and termed a "hedgehog" black hole. This was then 

■^Since this action actually describes a one-dimensional (radial) system, it may be possible to 
obtain a superpotential for it to explain the perfect square structure, and to derive first order 
equations. I thank A. Faedo for drawing attention to this point. 

^In pl[ , the backreaction of a uniformly smeared configuration of strings connecting antipodal 
points of the boundary was analyzed (in global AdS^). 
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generalized to asymptotically AdS spacetimes in |31]. The characteristic feature 
of the so-called hedgehog configurations in AdS^ is that the gravitational potential 
includes a term ^ —Q/r for large radial coordinate r, directly following from the 
linear dependence of the mass of a stretched string on its length. 

In the present case, in the gauge grr = jr, we obtain the following large r 
asymptotic expansions for the metric functions and the dilaton, 

. _ 1 . _ 2 Sg £ , 39233 

9rr-y,. ^« - ^ " 97";^+ 158760 7^ ' ^^-^^ 

,2. ^2 , 4g . 877 
9r 3r2 158760 

3r3 H 504 r6 

2 Q 5441 Q\ v^-e Q 
e^'? = 1 + — — ■ - — — — — H — H z + ... 



35 r3 98000 22 



We note that the leading correction to the AdS^ metric scales as ~ —Q/r, which is 
due to the macroscopic string source in the bulk. This breaks conformal invariance 
and Lorentz invariance of the boundary theory and the metric is not strictly asymp- 
totically AdS. However it is of the asymptotically locally AdS (AlAdS) form ||33| i.e. 
for large r, 

R^..x. = (g^agux - g^xgu.) {l + O (^)) . (3.2) 

After fixing the asymptotic form of the metric, the solution space is characterized 
by three integration constants e, f4 and fg corresponding to normalizable modes. 
The active fields in the solution are the dilaton 0, the metric component gu, and 
the s-wave components of the volume scalars a, i] of the and respectively. 
The volume scalar rj is dual to the irrelevant, dimension eight operator Os = 
-^STr [F"^ — |(F2)2], where 'STr' denotes the symmetrized trace |Q On the 
other hand the dilaton is dual to the marginal operator 04 = -^TrF^. The inte- 
gration constants V4 and vs are the expectation values of O4 and Os respectively, 
whilst e is related to the energy density of the boundary field theory. Note that we 
have taken the value of the dilaton (f) to vanish asymptotically. In principle this is 
a tunable (dimensionless) parameter in the solution corresponding to the marginal 
gauge coupling of A/" = 4 theory. 

In summary, the heavy quark density in the boundary gauge theory does not 
introduce irrelevant operators/non- normalizable modes in the bulk. It does mildly 
alter the asymptotics. In the absence of dimensionful VEVs V4, vg and e, the back- 
ground can be obtained as an expansion in powers of ^ since Q is a number density 
and has canonical scaling dimension three (in the UV). In general, we have no reason 
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to expect the VEVs to vanish. The existence of a smooth infrared (IR) solution 
may well require the dimensionful VEVs to be turned on (both at zero and non-zero 
temperatures). The possible values of these can only be ascertained after the IR 
dynamics of the finite density state and the complete flow towards it is understood. 

Before turning to the putative IR behaviour of the model under consideration, 
we briefly touch upon features of planar AdS-hedgehog geometries in Einstein gravity 
and how holographic regularization of the action of so-called hedgehog configurations 
proceeds. 



3.1 Holographic regularization 

A potential subtlety associated with the background above concerns the definition of 
a properly renormalized bulk action. This is an issue since the background contains 
sources which generate corrections (~ —Q/r) to the bulk AdS^ metric. To clarify 
this, let us first look at the divergent terms in the full 5D bulk (Euclidean) action. 



obtained by plugging in the IIB hedgehog asymptotics ( pyif) 



^1 



bulk 



I d'x I rfr(-8r3-f Q + ...), (3.3) 

where A is the UV cut off, which we treat as the boundary. Following the standard 
approach towards holographic renormalization |38|, we must also include the 

Gibbons-Hawking (GH) boundary term and possible geometrical counterterms. In 
the situation with a planar boundary the only non- vanishing counterterm arises from 
the boundary cosmological constant. So we obtain, 

Sgh = d'xVhTrKl^^ = -^Jd'x ( 8 A^ + | Q A + . . .) (3.4) 

?ct = — j d'x^Ql^^ = —J d'x (6A^ + |gA + ...) . 



87r2 

Here hap and K are the induced metric and the extrinsic curvature of the boundary 
at r = A, respectively. Summing the bulk and boundary actions yields 

^buik + SGH + S,t = ^ d^xAQA + ... (3.5) 

The usual boundary terms are sufficient to cancel the leading divergences, leaving 
behind a linear divergence proportional to the number density of strings/quarks. We 
will now see that this term can also be removed by adding an appropriate counterterm 
for the source action, namely a boundary term for the Nambu-Goto action of the 
macroscopic strings that source the geometry. In particular, a string stretching to 
the boundary of AdS space has to satisfy a Neumann boundary condition along 
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the radial direction, and four Dirichlet conditions along the directions parallel to 
the boundary. To this end, we need to include a term in the F-string action which 



replaces the radial coordinate with its conjugate momentum 



Spi Spi + SSfi = J dtr^^\^^^ . (3.6) 

Here r' denotes a derivative with respect to the world-sheet spatial coordinate. For 
the smeared configuration of strings this is precisely given by 

^Sfi = J rf^x4QAe^/2-io^/3y^|^^^, (3.7) 

which cancels the linear divergence of the bulk action. This divergence^ has a simple 
physical interpretation - it corresponds to the bare mass of the heavy quark, which is 
formally infinite and must naturally be subtracted to define a sensible energy density. 

3.2 AdS-hedgehog black hole in pure gravity 

In pure gravity in five dimensions with a negative cosmological constant, the back- 
reaction of a uniform distribution of strings gives rise to a background with metric 
(in Schwarzschild-like coordinates). 



ds' = -fif)dt' + ^^+f'dx\ /(r) = r-^ - ig - (3 
j[r) 3r r"^ 



This has a horizon at r = r^j. In a Fefferman-Graham type expansion (where the 
boundary is at r — )■ oo), the metric can be recast as 



-A-Qrh 



which agrees with the asymptotic behaviour eq.( |3.1| ). The thermodynamics of such 
solutions is non-trivial, particularly when the horizon is spherical In the planar 



case, it is easily seen that the hedgehog black hole in GR has a smooth zero temper- 
ature limit wherein /(f) has a double zero at f = r/i, resulting in a non- vanishing 
entropy density s. 



Q _ Q 



1 

^ 2.3- '"^^if- <'-'°' 

Thus the zero temperature system (in Einstein gravity) has a non-zero, large entropy 
density and the IR physics is described by the near horizon AdS2 x geometry. 
This feature, similar to the Reissner-Nordstrom black hole, seems unwanted and we 
may expect the type IIB setup dual to A/" = 4 SYM with heavy quark density to 
have a different and physically more appealing IR behaviour. 



^The prefactor fp: x AQ, in eq.(|j) is equal to 
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4. Exact (IR) scaling solution 



The most interesting question about the system we are studying is with respect to 
its low energy, long wavelength behaviour at low/zero temperature. We have seen 
that, in the pure gravity case, the low energy physics is dictated by a near horizon 
AdS2 region. We would now like to understand whether the type IIB hedgehog has 
a substantially different IR behaviour leading to a zero entropy ground state at zero 
temperature. 

We find that the system of "IIB plus strings" also admits an exact scaling solu- 
tion. It was not a priori clear that the system described by eq.{T^^) should possess 



such a solution. In fact, one may have naively expected the UV AdS-hedgehog ge- 
ometry to flow in the IR to a background sourced by multiple D5-branes wrapped 
on (corresponding to baryon vertices located at the origin of the space). How- 
ever, closer inspection reveals no clear separation of scales (in terms of the radial 
coordinate r) wherein the energy density in the three-form flux ~ dominates over 
the cosmological constant and the string density. This leads us to look for scaling 
solutions where all terms in the potential are comparable. We find, using the same 
coordinate system as in the UV hedgehog asymptotics, namely with cjrr = ^/r"^, an 
exact scaling solution to the equations of motion, 

dsl) = ^ - r^'/Ut' + r'IUx\ i = VW (^^y , (4.1) 

''''' e^^^ri^V'. 1.0297. 



578 734 Q2 ' ^121, 

This is a constant curvature spacetime, which we can cast in a more conventional 
form by defining a new radial coordinate p = r^/^, and after a rescaling 

dsl^ = + p^dx^ - p'uA , (4.2) 



2 144000 f, 
P 



17 11 Q2 

Therefore the scaling solution describes a regime with Lifshitz-like anisotropic scaling 



2^] and dynamical critical exponent z = 7. The scaling symmetry is not exact due 
to the logarithmic running of the dilaton with the radial coordinate p. 

The first observation is that the dynamical exponent is independent of the heavy 
quark density Q. It is also clear that the full ten dimensional metric is smooth in 
Einstein frame, since the has a constant radius e'^ = -y/ y|| ^ 1.015, while the 
Ricci scalar for the Lifshitz metric with z = 7 is, 

R = -2{z^ + 3z + 6) r 2 ~ -14.06 . (4.3) 
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However, as the coupling e'^ runs to zero at p = 0, the string frame metric has a 
curvature singularity in the deep infrared and the solution cannot be trusted since 
it will receive large a' corrections ^. Nevertheless, as we will see, we can always 
introduce a temperature into the system and cloak the singularity behind a horizon, 
without spoiling the scaling regime and therefore steer clear of the potential problems 
near p = 0. In this respect the solution ( ^.21 ) is quite similar to the charged dilaton 
black brane metrics of P21 E31 . 



Since the dilaton grows without bound for large r, in order to make sense of this 
solution we must be able to embed the scaling regime within the flow originating 
from the asymptotically (locally) AdS hedgehog background. 

It is worth noting that any putative baryon vertices (wrapped D5-branes) which 
could be regarded as sources of the F3 flux, are absent, possibly hidden by the 
"Lifshitz horizon" . Hence, for any small non-zero temperature, the system of heavy 
quarks is in a deconfined phase. Within the type IIB framework we can act on our 
system of smeared fundamental strings with the SL(2,Z) duality group. The action 
of S-duality, for example, replaces the F-strings with D-strings (corresponding to 
magnetic sources on the boundary) and simply changes the sign of the dilaton, so 
that e'f' diverges at p = 0, requiring the inclusion of quantum (loop) corrections in 
the deep IR. 



4.1 Effective Potential 

Let us now briefly explain the appearance of the Lifshitz solution by examining the 
equations of motion. The UV AdS-hedgehog solution nominally depends on two 
parameters: Q and the asymptotic value of the dilaton. However, in the absence 
of a scale in A^ = 4 SYM, bulk gravity solutions with different values of Q are 
related by a straightforward rescaling of the radial coordinate, and therefore do 
not really represent new solutions. On the other hand, the potential in the 5D 
action ( p.7| ) depends on the combination Q'^e^, so that UV solutions with different 
asymptotic values of the dilaton are related by a rescaling of Q. Given this, we may 
distinguish different UV solutions by the asymptotic value for the parameter e"^ . 
Independently of this asymptotic value, in the scaling Lifshitz-like solution eq.([4.2|), 
the combination e*^ is completely fixed and there are no free parameters (apart 
from coordinate rescalings). This indicates an attractor behaviour along the lines of 



what was found in [23, E4l. 



It is perhaps easiest to analyze the equations of motion in the gauge gu = V^'ee 
where q is an appropriate radial coordinate. From eq. (|2.7|), one identifies the natural 



^Note that this string frame curvature singularity is completely distinct from the effect of diver- 
gent tidal forces of backgrounds with exact Lifshitz scaling pq, Hl[ 
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effective potential for the one-dimensional (radial) problem, 

V = ^ {Q e<^/2 eioV3 ^-sa ^ 4 ^-20,73^2 _ ^3^-16,73 _ (4 4) 

The equations of motion (one of which is a constraint) in this gauge are, 

[gtt (e'^)')' + V = 0, {gu e=^^ r]')' ~ §^ d,V = , (4.5) 

The effective potential has a runaway behaviour, with a critical point at infinity 
where the dilaton e'^ and the volume scalar e"^ vanish (and so does the potential 
itself). The equations of motion then admit a scaling solution where this critical 
point is reached at the horizon {g = 0) by logarithmic running of (j) and a. In this 
scaling regime the scalars r] and the combination e'^~^'^, both remain fixed at constant 
values. It is a straightforward excercise to check in this gauge that the field equations 
are solved by 

72 , — 

n2 ^ »/ 

121 ' 



9u = ^q\ e^=\7Tif' (4-6) 



pS- _ 177 A _ 6 11979 677 



where i is defined in eq. ([4.1|) The individual normalizations of e'^ and e*^ can be 



changed by a coordinate rescaling subject to the requirement that the combination 
g0-6(7 fixed. After a coordinate transformation the background can be put in the 
Lifshitz-like form 



4.2 Black hole in scaling regime 



We have already pointed out that the dilaton running to zero in the infrared implies 
large string frame curvatures in the deep IR, requiring the inclusion of a' corrections 
in the gravity dual. Introducing a finite temperature would shield the system from 
this regime where the gravity description becomes inaccurate, and provide a con- 
sistent gravity dual of the boundary field theory at strong coupling. We find that 
the field equations admit an exact solution which is a simple deformation of the Lif- 
shitz background, yielding the following metric with a horizon in Schwarzschild-like 
coordinates: 



-gtt dt^ + ^ + ^1'^ dx"^ 
gtt 



(4.7) 



9tt = Q 



1-1^ 
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while the dilaton and the 5*^ radius remain unchanged. We therefore deduce that 
the Hawking temperature T and entropy density s for the black brane scale as 

r = 5|, .~r'/'(3»/'. (4.8) 

The dependence of the entropy density on Q is fixed by dimensional analysis. The 
scaling of the entropy with temperature is also exactly as expected for a theory with 
Lifshitz scaling symmetry {s ~ T^'^~^^^^ in d spacetime dimensions). Taken at face 
value this implies a vanishing entropy at zero temperature, although for vanishingly 
small T the classical solution should not be trusted. At any non-zero temperature 
the dilaton runs down to a non zero value oc T^/'' at the horizon. 



5. Flow from AdS to Lifshitz 

Finally, we would like to obtain evidence indicating that the Lifshitz-like scaling 
solution actually emerges as the infrared description of the N' = 4 theory with heavy 
quark density. To show this we would need to construct a flow connecting the AdS^ 
(hedgehog) background to the Lifshitz scaling regime. There is a natural dimensionful 
scale r ~ Q^^^, above which the geometry should approach the hedgehog geometry. 
Below this scale, all terms in the effective potential for the system compete with each 
other and the system should then enter the Lifshitz-like scaling regime. 

The natural way of constructing such a flow first requires the identification of 
relevant and irrelevant deformations of the putative IR background. We have shown 



that the UV AdS asymptotics ( pTip does not contain sources for any relevant de- 
formations (preserving the requisite symmetries), although it does allow VEVs for 
certain local operators. Perturbing the IR solution by an irrelevant deformation and 
applying a numerical 'shooting' method should enable one to reach the AdS^ asymp- 
totics in the UV. We parametrize the fluctuations 6i{g) about the Lifshitz solution 
as, 

9u = g^l + S,{g)) , a = '^Ing + S^ig), (5.1) 

At the linearized level, the fluctuations satisfying the equations of motion can each 
be expressed as a power series, with the leading terms given by, 

6i{g) ~ flip'', ^2 — CL2 +0.30"', ^3 — 602 + 04^"', ^4 ~ 05 ^".(5. 2) 

We find six possible values for n, 



-1 , n±± = ± i (917 ± 8v^) , (5.3) 



= ^, ri2 
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where all four sign combinations are allowed. Depending on the values of n, the 
coefficients in the series expansion satisfy special relations and can be identified as 
relevant and irrelevant perturbations. For each allowed value of n, there is precisely 
one independent coefficient. The perturbation with exactly the one that 

leads to the black brane metric and is relevant in the IR. The n = — 1 fluctuation 
corresponds to an additive shift of the radial coordinate. The two perturbations 
that are irrelevant in the IR and can be used to deform the Lifshitz background 
to alter the UV asymptotics are the ones corresponding to n = n++ = 1.973 and 
n = n+_ = 1.23. Using these perturbations, and adjusting their coefficients suitably, 
we have numerically found a flow interpolating between the Lifshitz scaling regime 
and AdS^ x 5*^ in the ultraviolet (see Figs. (|I|) and (^). 




- ln(r) 



Figure 1: Left: The 5D Ricci scalar changes from 152/^2 ~ -14.06 in the IR Lifshitz 
background to —20 in the UV, the latter corresponding to AdS^ asymptotics. Right: The 
radius of the 5^ flows from unity for large r to y ~ 1.014 in the IR. 



(p cr 




Figure 2: Working in the gauge Qrr = 1/?"^, the dilaton e'^ ~ ^6/£^ ^j^g volume scalar 
e'^ ~ r^/^ and the metric gu ~ r^^/^ in the scaling regime. Left: The dilaton runs according 
to the r^/^ power law in the IR and vanishes in the UV. Right: obeys the r^/^ law in 
the IR and is linear in r, as expected in the UV AdS region. 

The numerical solutions depicted in flgs. (P and (0) were obtained in the gauge 
g„ = ^/t"^ in which the UV and IR asymptotics take the forms shown in ( |3.1| ) and 
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( [4.1|) , respectively. The coordinates r and g are related as r = g^^'^. All plots were 
obtained for Q = 1.89 and — )• as r — oo. 

In principle, given the numerical flow it should be possible to read off the VEVs, 
V4 and Vs, of the dimension four and dimension eight operators, (^4 and Og respec- 
tively from the UV asymptotics ( |3.1| ). However, we have been unable to accurately 
isolate these higher order (~ and r~^) terms numerically. There is no reason 
to expect them to vanish. The fact that there are precisely two independent irrele- 
vant deformations of the IR Lifshitz point is consistent with the presence of the two 
VEVs in the UV solution. It should be possible to find a relation between the IR 
deformation parameters and the values of the V4 and vs ^. 



6. Discussion 



The main focus of this work was to understand the influence of a finite density of 
heavy quarks on the dynamics of a known gauge theory (namely JV = 4 SYM) with a 
strong coupling string dual. The fact that the system exhibits a flow to a Lifshitz-like 
scaling regime appears interesting and there are certain aspects of this that deserve 
further study: 



Since the microscopic description of the dual field theory is available, it po- 
tentially provides an opportunity to understand, in field theoretic terms, the 
emergence of the IR scaling regime and the breaking of it by the running dila- 
ton. The insertion of any number of straight timelike Wilson lines, induces a 
shift in the Lagrangian density, 

6C r^Yl ^'(^ - ^'^) - ' ^o(^' ^) - ' ^)) , (6.1) 

k 

where 9^^^-^ is a unit vector in fixing the coupling of the /c-th quark "fiavour" 
to the six scalars in A/" = 4 SYM (equivalently, the position of a macroscopic 
string on the 5*^). The Q(k) can be regarded as one-dimensional fermions living 
aX X = Xk (see e.g. Utilizing this as the starting point, it might be 

possible to shed light on the long-wavelength physics or the ground state of the 



system, by employing a Hartree type approximation, as in |^2[, encoding the 
collective effect of the heavy quark impurities. 

Perhaps the most natural question following from the gravity picture is whether 
the IR scaling behaviour is in any sense generic, i.e. can also be shown to arise 
from similar finite density configurations in other field theories with known 



thank A. Faedo for remarks on this point. 



-14- 



string/gravity duals (e.g. M = 2 and M = 1 SCFTs [43] in four dimensions 
with AdS^ X duals and the ABJM model in three dimensions [Q)- ^ related 
issue is how the value of the scaling exponent z depends on the microscopic 
field theory in question, or equivalently, on the corresponding active modes in 
the gravity dual. 

The smeared configuration considered in this paper preserves the 5*0(6) in- 
ternal global symmetry, but possibly breaks supersymmetry It should be 
possible to allow the strings to be localized on (breaking 50(6) to 5*0(5)) 
while smearing along the spatial directions, and subsequently obtain backre- 
acted geometries with some unbroken supersymmetry, similarly to p9l BO]. 



In the weakly coupled finite N D-brane picture, each individual heavy quark in 
our setup can be viewed as an Fl-string stretching between D3 and D5-branes, 
with 8 Dirichlet-Neumann boundary conditions. As is well known, such a 
string behaves as a fermion ^ . It remains to be seen if this fact makes 
the backreacted smeared string setup relevant for fermion physics in general, 
at finite density and strong coupling. 

In the context of the Af = 4 theory, the system with smeared heavy quark 
sources can also be formulated on 5^ [|l^. It is well known that A/" = 4 
SYM on 5^ exhibits non-trivial large- thermodynamics, both at strong and 
weak coupling p, The Hawking-Page deconfinement transition at infi- 

nite 't Hooft coupling has its zero coupling analogue which is driven by the 
Hagedorn growth of states. For theories with fundamental flavours on 5*^, in 
the Veneziano large- A^ limit, continuous deconfinement transitions have been 
shown to exist at finite quark density, and these are driven by the exponen- 
tially large degeneracy (~ e^) of baryon-hke states It would be extremely 



interesting to understand if there is any sign of this physics within the smeared 
string setup in global ^4^55. The hedgehog configurations in global AdS^ in 
pure gravity, are already known to exhibit a very interesting phase structure in 
the T — Q plane: a line of first order Hawking-Page transitions ending at a crit- 
ical point [^. Furthermore, for all values of Q, there exist zero temperature 
black holes with near horizion AdS2 region. 

Another motivation for investigating the above setup on 5*^ is provided by the 
IR cutoff given by the radius L of the three-sphere. In the range L Q^^^ ^ 1 we 
may expect the theory to approach the IR Lifshitz-like scaling behaviour, but 
eventually depart from it due to the finite infrared cutoff. This could be one 
way to regulate the deep IR behaviour of the scaling solutions with running 
dilaton. 



^We have not checked this, but it appears hkely since the individual strings/ Wilson lines have 
different orientations along the internal directions. 
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It is not difficult to see that the configuration studied by Headrick in [31] also 
admits a Lifshitz-like scaling solution in the planar case. The goal of that work 
was to obtain the free energy of A/" = 4 SYM as a function of the Polyakov 
loop order parameter i.e. an effective potential for the Polyakov loop at strong 
coupling. It would be interesting to understand what role, if any, is played 
by the low temperature scaling regime in obtaining a complete picture of this 
effective potential at strong coupling. 

A key missing ingredient in the heavy quark limit is a U{1)b global symmetry 
corresponding to conserved baryon number. As pointed out earlier, this global 
symmetry and its associated gauge field in the gravity dual will make an ap- 
pearance once the system is embedded within the framework of A^ = 4 SYM 
coupled to a large number of smeared fiavour D7-branes giving rise to mas- 
sive multiplets. For massive fiavours, the backreacted geometry in the smeared 
D7-brane setup of |^ (see also 0, |T^) has an IR geometry (below the quark 



mass scale) which is simply AdS^ x 5*^. Switching on an electric field on the 
D7-branes (corresponding to a chemical potential/quark density) will result in 
a spike on the D7-branes behaving precisely like a bundle of smeared strings in 
AdS^ X S^. Given the scaling behaviour seen above, it appears likely that em- 
bedding the smeared string hedgehog configurations within the D3-D7 setup, 
will lead to the same IR physics. The presence of a conserved U{1)b current 
would make it possible to define and compute transport properties of the IR 
theory. 
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